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A SCALING TECHNIQUE FOR AIR-BEARING SIMULATION 
OF PRECISION POINTING SYSTEMS 

By Gordon F. Bullock and Frederick R. Morrell 
Langley Research Center 

SUMMARY 

A scaling technique which can be used to simplify the laboratory simulation of the 
precision pointing system for a large space telescope is presented for the single -axis 
case. It is shown that the dynamic response of the scaled and unsealed control systems 
can be maintained the same by suitable scaling of the vehicle inertia, disturbance torque, 
and the attitude rate and position loops. The control system for the telescope consists 
of a single -axis, twin-rotor control moment gyroscope as the momentum exchange device, 
and a high-gain position loop including a star sensor. Experimental results from an air- 
bearing simulation of the pitch-axis precision pointing system for three different scale 
factors demonstrate the validity of the scaling technique. 

INTRODUCTION 

To evaluate the precision pointing capability of a large space telescope, a vehicle - 
control simulation was designed to study stability, error, and response of the system. 
Because of very large moments of inertia about its axes and the small environmental dis- 
turbance torques exerted upon the telescope, scaling techniques other than conventional 
are required for a laboratory simulation. Using conventional scaling techniques would 
require the disturbance torque to body inertia ratio to be maintained. In the large space 
telescope this requirement would cause the scaled disturbance torques to be in a range 
where simulator disturbances due to unbalance, turbine torques, and air currents in the 
room would cause excessive errors. The simulator would have to be stabilized to 
approximately 48.5 nrad (0.01 arc sec), which is beyond the state of the art in simulation 
today. For this type of problem, an error -scaling technique is described which results 
in a change on the torque scaling and on the disturbance torque to body inertia ratio. The 
features of this scaling technique are (1) the dynamic frequency response of the scaled 
system is the same as that of the unsealed system and (2) the inputs to the control 
actuator are unsealed and allow the use of full-scale hardware in the simulation. 



SYMBOLS 


D gyro damping, N-m/(rad/sec) (ft-lb/(rad/sec)) 

H angular momentum of control moment gyro rotor, N-m-sec (ft -lb-sec) 

Ig control moment gyro gimbal inertia, kg-m^ (slug-ft^) 

I m inertia as seen by control moment gyro servomotor, kg-m^ (slug-ft^) 

Iv telescope pitch axis moment of inertia, kg-m^ (slug-ft^) 

K s sensor gain, volts/rad 

k inertial scale factor, nondimensional 

N gear ratio 

s Laplacian operator, sec“l 

T c control torque, N-m (ft -lb) 

Id disturbance torque, N-m (ft-lb) 

T m torque motor output, N-m 

a control moment gyro gimbal angle, rad 

0 V telescope pitch axis attitude, rad 

t s sensor time constant, sec 

Dots over symbols indicate derivatives with respect to time. 

VEHICLE DESCRIPTION 

NASA Langley Research Center has examined some of the requirements associated 
with the design and operation of a large space telescope (refs. 1 to 4), an artist's con- 
ception of which is shown in figure 1. The telescope under investigation is a 3 -meter 
(120-inch) Cassegrain type that uses control moment gyros (CMG) as control actuators 
and a reaction control jet system for desaturation. In a 402-km (250 -statute -mile) 
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circular orbit the maximum disturbance torque exerted on the telescope is due primarily 
to gravity gradient. For the telescope configuration chosen, this torque is approxi- 
mately 0.298 N-m (0.22 ft -lb). The control system in its most stringent requirement 
must stabilize the telescope to better than 48.5 nrad (0.01 arc sec). The vehicle 
moments of inertia about its pitch and yaw axes are approximately 180 000 kg -m2 
(1.33 x 1 0 5 slug-ft2) and approximately 28 000 kg-m2 (2.07 x 10 4 slug-ft2) about the 
roll axis. The representative attitude control system investigated uses twin-rotor, 
single -degree -of -freedom control moment gyros as torque sources, with each rotor of 
the pair having an angular momentum H of 407 N-m-sec (300 ft -lb-sec), a gimbal 
inertia I g of 0.203 kg-m2 (0.15 slug -ft 2 ) and damping D of 0.854 N-m/(rad/sec) 

(0.63 ft -lb/(r ad/sec)). 

CONTROL SYSTEM DESCRIPTION 

A simplified block diagram of the pitch-axis attitude control system of the telescope 
is shown in figure 2. The single-axis rigid-body equation of motion of the telescope is 
(ref. 5) 


I v $v = T d " 2Ha cos a 


(1) 


A disturbance torque T^ acting on the telescope inertia I v produces changes in the 
telescope rate 6 y and attitude 0 y . This telescope rate is detected by the CMG gimbals 
through the gyroscopic cross -coupling torque. The telescope attitude is determined by 
the star sensor which transmits an error signal to the CMG gimbals. The error signal 
controls a torque motor to drive the CMG gimbals; this signal changes the CMG momen- 
tum vector which generates a control torque T c opposite to the disturbance torque. 
Analog computer studies (refs. 3, 4, and 6) have shown that a high loop gain is necessary 
to stabilize the telescope to 48.5 nrad (0.01 arc sec). Furthermore, these studies showed 
that nonlinearities in the CMG result in undesirable limit cycle operation. To check the 
validity of the assumptions made in the computer studies, a single -axis air -bearing sim- 
ulation of the attitude control system was designed. By using conventional scaling tech- 
niques, the problems associated with the large telescope moments of inertia and the small 
disturbance torques were immediately apparent. If conventional scaling techniques were 
used, the disturbance torque to body inertia ratio would be maintained. A simulation of 
the pitch-axis attitude control system, a simulator inertia of 40.7 kg-m 2 (30 slug -ft 2 ) 
being assumed, would require the disturbance torque to be 


^(simulator) = T^(telescope) 


Iv 

Iv 


(simulator) 

(telescope) 


0.298 40,7 , = 67 juN-m 

1.81 X 10 5 


(0.0095 in-oz) 


(2) 
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Vehicle 



Figure 2.- Block diagram of single -axis telescope attitude control system. 

This value of disturbance is approaching the limit of what can be reasonably achieved on 
a good air -bearing simulator and represents the maximum allowable disturbance torque. 
Furthermore, pointing errors of less than 48.5 nrad (0.01 arc sec) must be detected. It 
is doubtful that either of these requirements can be met in a room subject to even small 
temperature gradients and air motion. Even if the other requirements could be met, the 
CMG unit is a pure torque device and does not lend itself readily to torque scaling with- 
out affecting its nonlinear characteristics. These nonlinearities are one of the principal 
reasons for performing the hardware simulations and should be retained. These consid- 
erations led to the conclusion that the attitude control system for a 3 -meter diffraction- 
limited telescope could not be readily simulated in the laboratory to 48.5 nrad (0.01 arc 
sec) by using conventional scaling techniques. 

Scaling Technique 

A simplified block diagram shown in figure 3 describes an error -scaling technique 
which alleviates these problems and allows the attitude control system to be simulated. 
The block diagram is for the same basic control system as shown in figure 2, but, in this 
case, the disturbance torque to body inertia ratio has been changed. The vehicle moment 
of inertia I v has been decreased by a factor k and the disturbance torque Td has 
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Vehicle 



Figure 3.- Block diagram of scaled single -axis telescope attitude control system. 


been reduced by the factor \fk; thus, the telescope rate and position are scaled up by /k. 
The scaled single -axis vehicle equation of motion is 



The factor k is chosen so that the vehicle moment of inertia is reduced to values that 
can be easily implemented in the laboratory and yet keep the disturbance torque high with 
respect to noise sources such as simulator torques and air currents in the room. In this 
scaling technique it is important that the star sensor detects unsealed star image motions 
if actual hardware is to be used in the simulation. This requirement can be attained by 
changing the focal length of the star -sensor objective optics so that the size and motion 
of the star image on the focal plane would be the same as that in the unsealed telescope. 
This optical scaling is represented by the l/tfk gain reduction in the attitude sensor 
shown in figure 3, and results in a sensor output that is the same as that in the unsealed 
telescope. Thus the sensor gain becomes 


Sensor output = 


|/k(9 v V^= Sensor gain 


(4) 
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The angular momentum H of the gyro rotors has been decreased by the |/k. This 
decrease compensates for the increase in telescope rate 9 V since 



= H0 V cos a 


(5) 


which is the same as that in the unsealed telescope. This condition in conjunction with 
the sensor modification renders all the inputs to the CMG unsealed; therefore, the CMG 
gimbal rate and angle are identical to those in the full-sized telescope. The control 
torque T c from the CMG is reduced by l/k since the momentum H is reduced. This 
procedure yields the correct reaction torque to counter the scaled disturbance torque. 

The open-loop frequency response for the telescope control system is obtained from 
figure 2 as 


2HK S 

2H 2 s(-^-§ s 2 + s + l) (t s + 1) 




2H^ 


The closed-loop frequency response is 


( 6 ) 


0v 

T d 


d \d s + + 


2h2s S s2 + s + W + 1 ) + 2HK s 


(7) 


where the CMG gimbals have been assumed to be 0° and no compensation terms have 
been included in the expressions. Similarly, the scaled open-loop frequency response is 
found from figure 3 to be 


2H Ks 
l/k |/k 

2H 2 
k 

The scaled closed-loop frequency response is 


2H 2 /k/~ ' 1 2H 2 / k. 


s 2 + 


n v D/k' 


Is + 1 


(t s s + 1) 


l/k0 v 

Td 


s + 1 K t sS + !) 


2H 2 s 


GvI 


v x g/ 


V2H 2 / k, 


Is 2 + 


I V D/ k ' 
\2H 2 /k j 


s + 1 


. _ 2HK S 

(ToS + 1) + — # 

S ^k |/k 


( 8 ) 


(9) 


7 



It can be seen from equations (6) and (8) and equations (7) and (9) that the dynamic 
responses of the scaled and unsealed vehicle would be the same. 

In figure 3 it can be seen that the inputs to the control moment gyros are unsealed; 
this condition makes it possible to use control flight hardware in the simulation by suit- 
able adjustment of the CMG rotor speed. Since inputs to the attitude control system com- 
ponents are identical to those for the real case, the nonlinearities in these components 
will enter the simulation correctly. 

To obtain experimental data to verify the scaling technique, laboratory simulations 
were made at different moments of inertia or values of k and with a known gimbal fric- 
tion nonlinearity. If the scaling technique is valid, the limit cycle frequency for the dif- 
ferent runs should be the same. The telescope rate and attitude should be increased in 
magnitude by \[k whereas the gyro gimbal rate and torque motor output should remain 
identical in all runs. The equipment used in the simulation runs is described in the fol- 
lowing sections. 


Control Moment Gyro 

The control moment gyro, a photograph of which is shown in figure 4 and in block 
diagram form in figure 5, is the twin-rotor, single -degree -of -freedom type. This CMG 
was designed and built as a research tool with provisions made to vary the gimbal inertia 
and gyro damping. The effective gimbal inertia can be varied between 0.1357 and 
0.2714 kg-m^ (0.1 and 0.2 slug-ft^). Gyro damping can be varied between 0.61 and 
1.22 N-m/(rad/sec) (0.45 and 0.9 ft -lb /(rad/sec)). 

Each rotor has a momentum storage capability of 6.03 N-m-sec (4.45 ft -lb-sec) 
when driven by a 445-Hz power source. Provisions were made to change the momentum 
storage capability by varying the frequency of the drive power. The two rotors are 
mechanically coupled through a 1:1 gearing. A servomotor is used to drive the gimbals 
through an 80:1 gear head. Stall torque of the servomotor is 0.0102 N-m (1.45 in-oz), 
a flywheel of 0.001356 kg-m^ (0.001 slug-ft^) attached to the servomotor represents the 
inertia seen by the servomotor. A 0.085 N-m (12 in-oz) negator spring is used on one 
of the gimbals to take up the backlash in the 1:1 gearing and in the 80:1 gear head. 

Star Sensor 

The star sensor is the pyramid beam splitter type as shown in figure 6. Light from 
the star entering the lens is focused on the beam splitter. From the beam splitter the 
light is directed to the two photomultiplier tubes located on each side of the beam splitter. 
Angular motion of the star with respect to the sensor will cause more light to be directed 
to one of the photomultiplier tubes. The outputs of the two photomultiplier tubes are 
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Figure 4,- Control moment gyro unit. 


L- 71-532 


subtracted to get the error signal. The sensor was designed for an output of 2060 volts 
per radian, which is the sensor gain required by the model gyro rate loop. 


Vehicle Inertia Simulator 

The air bearing, which is a three -axis vehicle inertia simulator used in a single - 
axis configuration, is shown in figure 7 with the CMG and control -system hardware 
mounted on the payload table. The simulator is supported at the center by a ball-and- 
socket type of low-friction air bearing and the payload table is mounted on a shaft attached 
to the ball. A journal air bearing is used to restrict the motions of the simulator to those 
about its vertical axis. The moment of inertia about the simulator vertical axis repre- 
sents the pitch axis of the telescope in the simulation. Power and signals to the simu- 
lator are channeled through a mercury -bath type of slip-ring assembly. 
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Figure 5.- Block diagram of control moment gyro unit. 
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Figure 6.- Pyramid beam-splitter-type star sensor. 
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CONTROL SYSTEMS SIMULATIONS 

Three simulation runs were made with moments of inertia of 277 kg-m 2 
(204 slug-ft 2 ), 138.3 kg-m^ (102 slug-ft^), and 46 kg-m^ (34 slug-ff2). In the 277 and 
138.3 kg-m2 runs, moment arms were extended out from the payload table and mass was 
added to provide the desired moments of inertia. 

In the first run the moment of inertia of the simulator was set at 277 kg-m 2 
(204 slug -ft 2 ). The gain in the position loop was set to give a 9.7 jurad (2 arc sec) 
pointing accuracy with 0.225 N-m (0.166 ft -lb) of disturbance torque on the vehicle. 

The CMG was driven by 445 Hz power to develop a total momentum of 12.1 N-m-sec 
(8.9 ft -lb-sec). A mechanical friction was applied to the torque motor flywheel on the 
motor side of the 80:1 gear head, of 0.0056 N-m (0.8 in-oz). This mechanical friction 
was sufficient to cause the system to go into limit cycle operation. The results of a run 
under these conditions are shown in figure 8(a). These results show that the vehicle 
rate and position peak at ±97 prad/sec (±20 arc sec/sec) and ±97 jurad (±20 arc sec), 
respectively. The gyro gimbal rate peaks at ±0.0035 radian per second (0.2°/sec). The 
output of the torque motor reaches ±0.00846 N-m (1.2 in-oz). The frequency of the limit 
cycle is approximately 0.218 Hz. 

In the second run, shown in figure 8(b), the moment of inertia of the vehicle was 
reduced from 277 kg-m 2 (204 slug-ft^) by a factor of 2 to 138.3 kg-m2 (102 slug-ft 2 ). 

The star sensor gain was reduced by the ^2 instead of physically changing the sensor lens. 
The drive frequency for the spin motors was reduced by the \/2 from 445 Hz to 314 Hz 
lowering the total gyro momentum from 12.1 N-m-sec (8.9 ft-lb-sec) to 8.55 N-m-sec 
(6.3 ft-lb-sec). The gimbal friction level was kept at the same value. The results 
of the second run show that the telescope rate and position reached approximately 
±0.131 mrad/sec and 0.131 mrad (27 arc sec), respectively. These values compare 
well with the predicted values of 0. 137 mrad/sec and 0. 137 mrad, respectively. The 
gimbal rate and torque motor outputs remained approximately the same. 

In the third run the simulator moment of inertia was lowered by a factor of 6 from 
277 to 46 kg-m2 (204 to 34 slug-ft 2 ). The position loop gain was reduced by \/6 by 
changing the gain of the star sensor. The spin motor drive frequency was lowered by 
the \(6 from 445 Hz to 181.5 Hz, and thus the gyro momentum was reduced from 12.1 to 
4.93 N-m-sec (8.9 to 3.63 ft-lb-sec). The gimbal friction level was not changed. The 
results of this run are shown in figure 8(c). The telescope rate and position peak at 
approximately 0.228 mrad/sec and 0.228 mrad (47 arc sec), respectively. These values 
are approximately \/6 times the output for the 277 kg-m 2 case shown in figure 8(a). These 
results compare favorably with rate and position predicted values of 0.238 mrad/sec and 
0.238 mrad, respectively. The traces showing gimbal rate and output of the torque motor 
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